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Abstract 



A system of traps is considered, each containing a large number of Bose- 
condensed atoms. This ensemble of traps is subject to the action of an exter- 
nal modulating field generating nonequilibrium nonground-state condensates. 
When the frequency of the modulating field is in resonance with the transi- 
tion frequency between two different topological coherent modes, each trap 
becomes an analog of a finite-level resonant atom. Similarly to the case of 
atoms in an electromagnetic resonant field, one can create entanglement be- 
tween atomic traps subject to a common resonant modulating field generating 
higher coherent modes in each of the traps. A method is suggested for regu- 
lating entanglement production in such a system of multitrap and multimode 
Bose-Einstein condensates coupled through a common resonant modulating 
field. Several regimes of evolutional entanglement production, regulated by 
manipulating the external field, are illustrated by numerical calculations. 
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I. INTRODUCTION 



Entanglement is a notion that is important for quantum computation and quantum 
information processing [1-5]. Considering entanglement, one usually studies the systems of 
resonant atoms, possessing internal states, or spin assemblies [6-8]. As possible candidates 
for engineering entanglement, trapped ions have also been considered [1,2]. On the other 
hand, recent advances in atomic physics with neutral atoms, experiencing Bose-Einstein 
condensation [9-14] , have led to several suggestions of employing these condensed atoms for 
producing entanglement. Thus, entanglement was shown to be feasible in Bose-condensed 
systems of atoms with two different internal states or with two collective modes [15,16]. 
Entanglement can be realized between the atoms from two clouds with different momenta, 
which arise in the process of coherent scattering of light from a Bose-Einstein condensate 
[17-20]. It is also possible to consider entanglement in Fermi-Bose mixtures [21-23] and 
in atom-molecule mixtures created by Feshbach resonance [24,25]. A review of the general 
properties of two or more entangled modes can be found in Ref . [26] . 

A promising and efficient way of generating several topological coherent modes in Bose- 
condensed systems is by means of the resonant modulation of the external trapping potential 
[27-30]. The properties of these coherent modes have also been studied in a number of works 
[31-45]. A Bose-condensed system with so generated coherent modes possesses a rich be- 
haviour displaying several nontrivial effects, such as interference fringes [29,41], interference 
current [29,41], mode locking [27-29], dynamic transitions and critical phenomena [28,29,37], 
chaotic motion [30,44], harmonic generation and parametric conversion [30,44], and atomic 
squeezing [29,43]. It was mentioned [42,43] that a trap with several generated coherent 
modes is an analog of a finite-level resonant atom. 

In the present paper, we consider not just a single trap, but an ensemble of many traps 
subject to the action of a common resonant modulating field, generating topological coher- 
ent modes in each of the traps. Such a system is equivalent to an ensemble of many resonant 
atoms subject to a common electromagnetic field. As is known [6], the atoms, resonantly 
interacting with a common field, become entangled, which is termed resonant entanglement. 
In the same way, there appears entanglement between the traps subject to a common res- 
onant field generating coherent modes. This resonant entanglement is a particular case of 
entanglement between coherent modes [26]. One of possible ways of creating a multitrap 
Bose condensate could be by means of optical lattices, which should be sufficiently deep, so 
that a potential well around each lattice site would contain a large number of atoms, thus, 
representing a separate trap. Then the lattice as a whole should be shaken by an external 
modulating field, resonantly generating coherent modes in each of the lattice-site wells. In 
an experiment, the resonant field can be produced by a time-dependent modulation of mag- 
netic fields or by laser beams [27-31]. The advantage of employing a deep optical lattice as 
a multitrap device is the possibility of varying the barriers between lattice sites and, hence, 
the feasibility of tuning the properties of the overall system. 

A peculiar feature of entanglement, occurring in multitrap multimode Bose condensates, 
is that this is an entanglement between mesoscopic (or almost macroscopic) objects, as 
compared to the entanglement between atoms, which are microscopic objects. At the same 
time, a mesoscopic Bose condensate is a quantum object, described by a wave function 
satisfying a nonlinear Schrodinger equation. This unique combination of being mesoscopic 



2 



in size and quantum in its features makes Bose-Einstein condensate a specific matter. A 
multitrap collection of multimode condensates constitutes a novel system that can find 
applications in quantum computing and information processing. 

Before going to a mathematical description, let us present in more detail the physical 
picture we keep in mind. A set of Bose-Einstein condensates can be created by loading 
ultracold atoms into a deep optical lattice. Thus a one-dimensional array of ^^Rb condensates 
was formed [46], with a period of 2.7 x 10~^ cm. The lattice was sufficiently deep, with the 
depth V^) 27r X 50 kHz a; 3 x 10^ Hz, or, in units of the recoil energy, Vq/ pa 600. The 
total number of lattice sites was 30, each containing ~ 10^ atoms. A similar experiment 
for ^'''Rb condensates was reported in Ref. [47]. A deep one-dimensional optical lattice, 
with adjacent sites spaced by 5.3 x 10~^ cm, contained about 10^ atoms. The number of 
lattice sites could be varied between 5 to 35, so that the number of condensed atoms in 
each site was between 100 to 10^. In these experiments [46,47], the number of condensed 
atoms N in each of the lattice sites was very large, reaching N ~ 10^. Each lattice site 
played the role of a trap containing a large number of completely Bose-condensed atoms 
in a well defined coherent state. But, because of the great depth of each site, the intersite 
tunneling was strongly suppressed. Such a state of matter is commonly called the Mott 
insulator. A specific feature of the insulating phase, formed in the experiments [46,47], 
is an enormous filling factor, with the number of atoms per lattice site being very large, 
N ^ 1. These insulating systems were probed [46,47] by releasing the atoms from their 
lattice, so that the wave packets of different condensates, emanating from separate lattice 
sites, could expand and overlap. A surprising discovery of these observations was that, 
contrary to what could be naively expected, there existed high-contrast interference fringes, 
which is commonly associated with the presence of coherence between different lattice sites. 
The interference-fringe contrast was indeed comparable to the results observed in the early 
interference experiments for two well-correlated condensates (sec discussion in the review 
papers [10,14]). Similar results were also obtained for deep optical lattices with low filling 
factors, showing that some kind of coherence persists even deep in the insulating phase, 
preserving a good visibility of the interference patterns observed after free expansion [48]. 
The theoretical explanation of these interference effects shows that, though separate lattice 
sites may look to a first glance independent, the overall insulating phase is a highly correlated 
state of the many-body system, being a generically entangled state [49]. 

An additional degree of entanglement for such deep lattices, with large filling factors, 
is provided by subjecting the whole system to the action of an external resonant field, 
being common for all lattice sites. Since the resonant field is common for the total lattice, 
one has to treat the insulator plus the field as a whole, as a single complex object. This 
situation is completely analogous to what one has subjecting an ensemble of finite-level 
atoms to a common resonant field, realizing the so-called resonant entanglement [6]. There 
is a straightforward and direct analogy between these two cases [41,42,50]. 

A Bose-condensed gas can also be made equivalent to a finite-level system. For this 
purpose an external alternating field, acting on the Bose-Einstein condensate of trapped 
atoms, has to oscillate with a frequency uj close to the transition frequency uoq between the 
ground-state level and that of the chosen excited mode. Because of the nonlinearity of the 
stationary Gross-Pitaevskii equation, the spectrum of self-consistent collective modes is not 
equidistant. This spectrum, for different trap shapes, was calculated and analysed in Refs. 
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[10,27,34,38,39,42]. Applying several driving fields, with the appropriately tuned frequencies, 
a multimodc system can be formed. Such a finite-level condensate is completely analogous 
to a finite-level atom subject to the action of a resonant electromagnetic field [41,42]. In 
the process of the resonant excitation of a Bose condensate, similar to the same procedure 
for an atom, there exist as well nonresonant transitions. However, their probability is much 
lower than that of the chosen resonant transition. The effect of nonresonant transitions 
was studied in Ref. [29], where it was shown that, for a typical setup, the occurrence of 
nonresonant transitions would destroy the condensate during the time approximately equal 
to the lifetime of atoms in a trap. Therefore, the resonant generation of nonground-state 
modes in a Bose condensate is a feasible procedure. These conclusions were obtained by 
employing accurate analytical calculations as well as by direct numerical simulations for the 
temporal Gross-Pitaevskii equation [30,44,45]. The resonant generation of a dipole mode 
was realized in experiment [32]. Vortices are the examples of the topological modes, which 
can be generated in the resonant way [31]. 

Each lattice site of a deep optical lattice represents a trap that may contain a large 
number of Bose-condensed atoms. According to experiments [46,47], this number can be as 
large as ~ 10"^. A collection of the lattice sites composes a strongly correlated system. 
The correlation between different sites comes from three sources. First, the standard method 
of preparing the considered system is by initially condensing a large atomic cloud and then 
imposing on it a periodic lattice formed by laser beams. Hence, at the initial time all atoms 
compose a general condensed system in a single coherent state, after which they are loaded 
in different lattice sites [46-48]. Second, though the tunneling between deep lattice sites is 
small, but, nevertheless, it is finite, because of which phase coherence on short length scales 
persists even deep in the insulating phase [48]. Finally, when the overall lattice is subject to 
a common resonant field, the whole system of the lattice plus field must be treated as one 
complex object. 

We show below that not merely one is able to generate entanglement in a multitrap Bose 
condensate, but, moreover, it is possible to effectively regulate the amount of the prodiiccd 
entanglement by manipulating with the applied resonant field. The potentialities of reg- 
ulating the entanglement generation is what makes the considered method important for 
applications. Using numerical calculations, we present several regimes of regulated entan- 
glement production, which illustrate wide possibilities of the suggested approach. In the last 
Section V, we discuss how the entanglement in the system can manifest itself in experiment. 

Throughout the paper, the system of units is employed, where the Planck constant and 
Boltzmann constant are set to unity, ^ = 1, /c^ = 1- 

II. SINGLE-TRAP COHERENT STATES 

Before studying a collection of many traps, let us, first, give an accurate mathematical 
description of coherent states in a single trap. The intention of the present section is to 
introduce the definitions and notations characterizing a single trap, which will be used in 
the following section for a generalization to the case of many traps. Let a trap contain 
A^ Bose atoms at a very low temperature, when the whole system is Bose-condensed. The 
Bose-Einstein condensate, as is known [10,51], corresponds to a coherent state jry > of a Fock 
space. The coherent state is an eigenstate of a destruction field operator, whose eigenvalue is 
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given by a coherent field 77(r) . In the Fock space, a normahzed coherent state is represented 
by a column 



\r]> = 

with respect to the index n = 0, 1, 2, . . ., enumerating the column rows. Equation (1) is the 
standard notation for a state in the Fock space [51]. A more detailed explanation of the 
notations and definitions, used in the paper, are given in Appendix. The scalar product of 
two coherent states is defined as 

oo g-AT 

where the scalar product of two coherent fields is 

(?7i,?72) = j rjl{r)r]2{r) dr . (3) 
The normalization conditions for the coherent states and fields are 



-N/2 



i=l 



(1) 



<77|77> = 1, {rj,r])^N. (4) 

A coherent state is a ground state for a Bose-condensed system [10,51,52]. Here we 
keep in mind a dilute system of atoms confined by means of a trapping potential U (r) and 
interacting through the local potential 

$(r) = $o5(r) , $0 = 47r ^ , 

m 

in which Ug is a scattering length and m, atomic mass. The stationary coherent fields are 
defined by the Gross-Pitaevskii equation 

H[r]k] r]k{r) = Ek r]k{r) , (5) 

with the nonlinear Schrodinger Hamiltonian 

H[v]^-^ + Uir) + %\ri\\ (6) 
zm 

The eigenproblem (5) yields a spectrum of the energies E^, labelled by a multi-index k. The 

lowest energy Eq = min^ Ek corresponds to the usual Bose-Einstein condensate, with the 
related coherent field i]o{r) being just the condensate wave function. The higher energies 
Ek correspond to the nonground- state condensates, with the wave functions %(r) called the 
topological coherent modes [27-30]. Each mode is normalized as 

(%,%) = AT. (7) 

The modes are named topological, since they describe topologically different distributions 
of atoms |r7fc(r)p. Recall a vortex that is the most known example of a topological mode. 
For each coherent mode r)k{r), one may construct, similarly to Eq. (1), a coherent state 



5 



\Vk > = 

Equation (8), in the same way as Eq. (1), is the standard notation of a state in the Fock 
space. The right-hand side here denotes a column with an infinite number of rows enu- 
merated by an index labeUing the generic column element inside the square brackets [51]. 
These mode coherent states are normalized to one, < rik\i]k >= 1, by analogy with Eq. (4). 
The coherent modes of Eq. (5), being defined by a nonlinear Schrodinger equation, are not 
necessarily orthogonal. Therefore the coherent states (8), with a scalar product 

< VklVp > = exp{-Ar + {rjk, %)} , (9) 

look also nonorthogonal for k ^p. 

Prom the Cauchy-Schwartz inequality we know that 

i(^?fc,^?p)r < {vk,vk){vp,vp) , 

where the equality takes place if and only if r]k{r) = '?/p(r). The labelling of the coherent 
modes is assumed to be done so that to avoid the degeneracy in the notation. That is, rik{r) 
can be equal to ?7p(r) then and only then, when k coincides with p. Hence, for different k 
and p, one has 

\{r]k,r]p)\ < N {k^p). 
Therefore, in the limit of a large number of particles, for the scalar product (9), we have 

lim <rjk\rip> = 4p • (10) 

iV— >oo 

This is equivalent to saying that the states (8) are asymptotically orthogonal for N ^ 1. 
Recall that here is the number of atoms in a single lattice site, which as is explained in 
the Introduction, can be made in experiments very large. 

It is worth stressing that the coherent states (8), composed of the coherent fields ?7fe(r), 
characterize collective states of an A^-atomic system, and should not be confused with in- 
ternal states of an atom. To avoid such a confusion, we deal here with atoms possessing no 
internal states. 



-N/2 



(8) 



III. MULTITRAP COHERENT STATES 

Now let us consider an ensemble of L traps, each of which can possess M modes, 

^ = M = (11) 

j k 

The traps are enumerated with an index j = 1, 2, . . . , L. As is discussed in the Introduction, 
such an ensemble of traps could be realized as an optical lattice, with the potential wells 
at each lattice site deep enough to incorporate a large number of atoms A^. For sufficiently 
deep potential wells, the condensate wave functions are localized close to the related lattice 
sites [10,53], so that the coherent modes ?7ifc(r) and 77jfe(r) of different traps practically do 
not overlap. 
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(12) 



This corresponds to an insulating lattice phase. For each coherent mode riik{r), we construct 
a coherent state \r]ik >, by analogy with Eq. (8). These coherent states, because of Eqs. 
(10) and (12), are asymptotically orthogonal, such that 

< ViklVjp > - ^iAp (A^ » 1) • (13) 
For each set {\r]jk >} of the coherent states \r]jk > a resolution of unity 

J2hk><Vjk\ -i 

k 

is asymptotically {N S> 1) valid, understood in the weak sense [10,52]. Thus, the set {\r]jk >} 
forms a basis, possibly overcomplete. 
Let us define a Hilbert space 

n,^c{\v,k>} (14) 

as a closed liner envelope over the set {\f]jk >}• Then, Tij is a Hilbert space of the coherent 
states associated with a j-th trap. The total Hilbert space for an ensemble of L traps is the 
tensor product 

n^^Hj, (15) 

3 

in which j = 1,2, ... ,L. The total space (15) contains entangled as well as disentangled 
states. We define the set of disentangled states 

P ^ {(g) \ip, > : >e n,} (16) 
j 

as a subset of H, composed of all disentangled states, having the form of the tensor products 
<S)j \<Pj>. 

When we deal with an absolutely equilibrium system, then, as is evident, all atoms 
condense to the lowest-energy state. This implies, assuming that all traps arc identical, that 
in each of them solely the states \r]jo > will be occupied, all other excited states being empty. 
In order to physically realize a nontrivial situation, it is necessary to invoke a mechanism for 
the creation of several coherent modes. Then, clearly, we have to consider a nonequilibrium 
system. An efficient mechanism for generating coherent modes has been suggested [27-30], 
based on a resonant modulation of the trapping potential. Dealing with a nonequilibrium 
system, the coherent states \f]{t) > from the Hilbert space (15) become functions of time t, 
generally, having the form 

\v{t)>^ Ecfe(^)(8)l^ifc> ■ (17) 

k j 

Here we assume that the same resonant pumping field acts on all traps, generating the 
coherent modes with the probability |cfc(t)p. 
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Varying the coefficients Ck{t) in Eq. (17), different entangled states can be created. 
Maximal entanglement is achieved when = 1/M. Thus, for |cfcp = 1/2, L = 2 and 

M = 2, Eq. (17) would represent a Beh state. For |cfcp = 1/2, L > 2 and M = 2, this 
would be a multicat state. And for |cfcp = 1/M, L > 2, with M > 2, a multicat multimode 
state would arise. 

The statistical operator, corresponding to the state (17), is 

p{t) = \v{t)><v{t)\. (18) 

To quantify the entanglement generated in the system, we shall employ the measure of 
entanglement production introduced in Refs. [54,55]. The advantage of this measure is 
its generality, which allows for the usage of the measure for any multiparticle systems; for 
example, the four-particle entanglement [56,57] can be easily characterized. The general 
measure of entanglement production [54,55] can be defined for an arbitrary operator. Here 
we shall use it for the statistical operator (18). 

Let us define a nonentangling counterpart of the operator (18) as 

m^<S)m, (19) 

where the partite operators are 

Pj{t)=Tr^n.^.}p{t)- (20) 

Entanglement, generated by the statistical operator (18), is quantified by the measure of 
entanglement production 



^ ml ' '''' 

in which the logarithm is to the base 2 and 

wmwv^snvwmfw (ii/ii = i) 

fev 

is the operator spectral norm, associated with the corresponding vector norm, and restricted 

to the set of disentangled states (16). 

For what follows, we shall need the notation 

nk{t) = \ck{t)W (22) 

characterizing the mode probabilities, or the fractional mode populations. As the partite 
operators (20), we get 

Pj(t) = Z]"'fe(*)l%fc >< Vjk\ ■ 

k 

Also, we find the norms 

■D = \\Pj{t)\\v = supnfc(i) , \\p®{t)\\v = [supnfe(i)]^ . 

k k 
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Then the measure (21) of entanglement, generated by the statistical operator (18), writes 
as 

e{p{t)) = (1 - L) log sup rikit) . (23) 

k 

This measure, for the present case, has a nice property of being expressed by the equation 

aim) = {L- l)e2{t) (24) 

through the bipartite measure 

£2{t) = - log supnfc(t) . (25) 

k 

The relation (24) here is valid for any number of traps L and any number of modes M. 
This relation is a particular concrete example of the fundamental property, formulated by 
Linden, Popescu, and Wootters [58,59], that "the parts determine the whole in a generic 
pure quantum state". 



IV. REGULATING ENTANGLEMENT PRODUCTION 

Entanglement is such a general effect that it exists practically for all condensed-matter 
systems, changing together with phase transitions [54,55]. However, in the majority of cases, 
it is not easy to govern the level of entanglement. To be practically useful, the system should 
allow for an efficient way of manipulating with the generated entanglement, as is the case 
of some simple spin systems [8] . 

Here we show how entanglement production can be regulated in a multitrap Bose- 
condensed system. To excite the higher coherent modes, it is necessary to use the method 
of resonant generation [27-30]. For this purpose, the system is subject to the action of an 
external modulating field 

y(r, t) = Vi{r) cosa;^ + V2{r) sina;i , (26) 

whose frequency is tuned close to the transition frequency 

uJo = Ei- Eo (27) 

between the ground-state mode, with an energy Eq, and an excited coherent mode, with an 
energy denoted by Ei. The resonance condition 



Aa; 



CO 



< 1 (Aa; = uj-uJo) (28) 



is assumed to be valid. For a nonequilibrium system, one has to consider the time-dependent 
Gross-Pitaevskii equation 

i^^{r,t) = {H[^] + v) ^{r,t), (29) 
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in which the wave function is normahzed to one, = 1, the nonhnear Schrodinger 

Hamiltonian is 



H[p] = - ^ + + ^oiVl^^r , (30) 

and V — V{y, t) is the resonant field (26). We look for the solution of Eq. (29) in the form 

(^(r,t)^^c,(t)(^,(r)e-^^'=*, (31) 

k 

where the coefficients Cfc(t) change in time slower than the exponential exp{—iEkt). The 
quantity |cfc(t)p defines the probability of generating a k-th mode, or the fractional mode 
population. For these probabilities, the normalization condition holds 

k 

Since we have assumed that all traps are identical, and the same resonant field acts on each 
of them, the probabilities |cfc(t)p do not depend on the trap index. 

A resonant field of type (26) generates one higher coherent mode. Wishing to create sev- 
eral topological coherent modes, one should apply several resonant fields with the frequencies 
tuned close to the transition frequencies of the chosen modes [30,44]. Here we consider a 
two-mode case, one of them being the ground-state mode and another, an excited coherent 
mode. 

Substituting the form (31) into Eq. (29) and employing the averaging technique, we can 
derive the equations for the coefficients Ck{t). This derivation, with all details and with a full 

mathematical foundation has been given earlier [27-30]. Below, we present only the results 
of this procedure. We shall need the notation for the interaction transition amplitude 

akp = ^oN 2\ipp\^ - \ipk\^) (33) 

and the pumping-field transition amplitude 

Pkp = iVk: B(pp) , (34) 

where B = Vi(r) — iV2{r). For the two- mode case, it is convenient to introduce the average 
interaction amplitude 

a = ^ («oi + aio) (35) 
and to write the pumping-field amplitude in the form 

Poi^Pe'\ (3=\Poi\. (36) 
We denote the effective detuning as 

5 = A + i (aoi - Q!io) , (37) 
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where 



A = Aa; + aoo ~ ctii • (38) 
The fractional mode amphtude can be written as 

Ck{t) = \ck{t)\exp{i7rk{t)t} , (39) 

where nk{t) is a real function representating the mode phase. The final equations can be 
cast to the system of two real equations for the population difference 

s{t) = \c,{t)\' - \co{t)\' (40) 

and for the effective phase difference 

X{t) = TTiit) - TToit) + ^ + A . (41) 



These equations are 



ds 

— — —(5\/l — sinx , 
at 



— = as -\ — , cos X + . (42) 

at vl — 

By scaling the time variable with a, one can see that the behaviour of the solutions to Eqs. 
(42) is defined by two dimensionless parameters 

b^^, e^^. (43) 
a a 

The former parameter characterizes the strength of the pumping field, and the second pa- 
rameter describes an effective detuning from the resonance, which can always be made small, 
l^l <^ 1. By solving Eqs. (42), one finds the fractional mode populations 

1 - sit) , , 1 + s(t) 

no{t) = , n,{t) = . (44) 

Finally, the entanglement generation is quantified by the measure (25) which, for the con- 
sidered two-mode case, is 

S2(t) = -log2sup{no(t),ni(i)} . (45) 

This measure changes with time in the interval < e2{t) < 1, being maximal when uq = 
ni = 1/2. 

Solving Eqs. (42), we suppose that at the initial time t — solely the ground-state mode 
exists, so that s(0) = —1 and x{0) — 0. Hence, initially, no(0) = 1 and ni(0) = 0, because 
of which £2(0) = 0. Then the pumping resonant field is switched on, which generates an 
excited coherent mode and produces entanglement. There are two ways of governing the 
amount of entanglement production. 



11 



First, one can vary the amplitude and frequency of the resonant pumping field, choosing 
by this the required parameters (43) , whose values define two main regimes of an oscillatory 
behaviour of the population difference s{t), and so, the regimes of S2{t). These arc the 
mode-locked and mode-unlocked regimes [27-30,34,37,41-44]. We solve numerically Eqs. 
(42) and calculate the measure of entanglement production (45). Keeping in mind that the 
detuning can always be made small, we set S — 0. Then the value be — 0.497764 is the 
critical point for the change of the dynamical regimes. Below be, the measure (45) oscillates 
with time, never reaching one. When the parameter 6, defined in Eq. (43), equals the critical 
point b — be, then the oscillating e2{t) reaches one. For the dimensionless amplitude of the 
pumping field b > be, the oscillations of S2{t) are always in the interval between and 1. 
But the oscillation period sensitively depends on the value of b. Thus, for b = 0.5, just a 
little above be, the period of oscillations is more than twice shorter than that for b = be- 
The period for 6 = 0.7 is about eight times shorter than for b — be- Thus, by varying the 
amplitude of the pumping resonant field, we can strongly infiuence the evolution of S2{t) 
both in its amplitude and period of oscillations. 

There is one more very interesting way of regulating entanglement generation, which can 
be done by switching on and off the applied resonant field. Recall that this alternating field 
can be easily produced by modulating the magnetic field forming the trapping potential 
in magnetic traps or by varying the laser intensity in optical traps. Then, it is possible 
to create various sequences of pulses for e2{t). For example, by switching on and off the 
resonant field in a periodic manner, we may form equidistant pulses of e2{t), with all pulses 
having the same shape, as is demonstrated in Fig. 1. But we can also switch on and off the 
pumping field at different time intervals, thus, forming nonequidistant pulses, as is shown in 
Fig. 2. The possibility of creating very different pulses is illustrated in Fig. 3. Regulating 
entanglement production by means of a manipulation with the resonant pumping field, it is 
feasible to organize a kind of the Morse alphabet. 

V. RELATION TO TIME-OF-FLIGHT EXPERIMENTS 

It is important to add an idea how the entanglement generated in the system could 
manifest itself in order to verify it in a measurement. This can be done by observing 
the expansion of atomic clouds released from the optical lattice. The absorption images 
and interference fringes, observed in such time-of-flight experiments, are known to strongly 
depend on the initial conditions of the atoms before their release. If the initial spatial 
distribution of atoms in one case is essentially different from that in another case, the 
interference pictures for such two cases will be noticeably different. The generation of 
topological coherent modes, from one side, is directly related to the level of the induced 
entanglement and, from another side, drastically changes the spatial distribution of atoms. 
Recall that the considered modes are termed topological exactly because of their principally 
distinct from each other spatial shapes [27-30]. Radically differing initial conditions, due 
to differently prepared modes, will inevitably result in distinguishable one from another 
absorption images. As is seen from Eq. (17), if there exists just one mode, say the ground- 
state mode, so that no — \ck\'^ — 1, hence the fractional populations for all other modes are 
zero, then \ri(t) > in Eq. (17) is a product state, with no entanglement. Respectively, the 
measures (24), (25), and (45) are identically zero. By generating excited coherent modes. 
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one makes the state (17) entangled. For instance, in the case of two modes, the maximal 
entanglement is achieved for \ck\ = l/V^, when Eq. (17) represents a multicat state. Then 
the measure of entanglement production (24) reaches its maximal value (L — 1) log 2 = L — 1, 
assuming that the logarithm is to the base two. Then the measures (25) and (45) are also 
maximal, being equal to one. 

Suppose that at the time to atoms were released from their traps. In the time-of-flight 
experiments one measures the atomic distribution 



in which p — mr/ {t — to) and p(r, r', to) is the first-order density matrix of atoms just before 
their expansion. The distribution (46) is normalized to the total number of particles NL, 
with N being the number of atoms in each lattice site and L, the number of sites, so that 



where p(r, to) = p(r, r, to) is the density of atoms before their release. 

The field operator, acting on the total Hilbert space (15), can be represented as a direct 
sum 



of the operators ipj{r) acting on the spaces Tij. By the definition of the coherent states, one 
has 




(46) 




(47) 



(48) 



V'jWI^fe > ^ r)jk{r)\Vjk > ■ 



(49) 




(50) 



Prom here, using Eqs. (48) and (49), we find 



p(r,r',to) = E E \ckito)\\kir) vUr') 



(51) 



k ij 



Substituting this into Eq. (46), we get 



pip) = E E Mto) VikiP) Vjkip) , 



(52) 



k ij 



were the Fourier transform 




(53) 



is introduced. Defining also the partial distribution 



Pfc(p) = E vUp) Vjkip) , 



(54) 
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we finally obtain 



P(p) = Mto) Pk{p) ■ (55) 

k 

When there is a sole mode, say when no (to) = 1, then there is no entanglement, and the 
measured distribution of atoms in the time-of- flight experiment is Po(p)- The generation of 
excited modes creates entanglement and at the same time essentially changes the distribution 
of expanding atoms (55). 

To give a more explicit illustration, let us consider a one-dimensional array of lattice 
sites in the tight-binding approximation. The ground-state wave function of atoms in a j-th 
site can be described by the localized orbital 

in which aj = ja, with a being a lattice spacing and j = 0, ±1, ±2, . . .. The effective width 
of the atomic cloud, given by Iq, is generally defined through a variational procedure taking 
into account atomic interactions [10,27,42]. The Fourier transform of Eq. (56) is 

= {ATif^4^, exp [- ^ pHl - ipa^ . (57) 
For the corresponding distribution (54), we get 



Po(p) = 2V7r NLloexp (-p'^ll) ^1 + XJ cospaij , 



(58) 



where Ojj = ai — aj. This is what one would observe in the time-of- flight experiment, when 
all atoms before the release are in the ground state and there is no entanglement. 

Suppose now that, in addition to the ground-state mode, the first excited mode is gen- 
erated having in a j-th site the wave function 



Here li is an effective width, which can be defined through a variational procedure [10,27,42]. 
Calculating the Fourieir transform of Eq. (59), we obtain 

%i(p) = -i2T:^/^yjNhpli exp (- ^ pHl - ipa^ . (60) 
Consequently, the related atomic distribution is 

pi(p) = 40F NLpHl exp (-pHt) + I E cospa^,- j . (61) 

The total observable atomic distribution in the two-mode case is the sum 
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p{p) = no{to)po{p) + ni{to)pi{p) 



(62) 



If both no (to) as well as ni(to) are nonzero, then entanglement is generated in the system. 
At the same time, the measurable distribution (62) is essentially different from the situation 
when just one term po{p) is present. In the particular case of two sites, when L — 2, the 
terms in Eq. (62) are defined by 

Po{p) = 4\/7r Nlo exp (-p^/o) (1 + cos pa) , (63) 

where a is the intersite distance, and 

pi{p) = Sa/tF Np'^lfexp {-P^lf) (1 + cos pa) . (64) 

It is evident that the sum (62) of the distributions (63) and (64) is drastically different from 
the sole term (63). This difference is especially noticeable at the point p = 0, where 

p(0) = no(to)po(0) = 8V^ Nlouoito) . (65) 

Value (65) diminishes as soon as the second mode is generated, so that no (to) becomes less 
than one. 

In conclusion, a multitrap ensemble of multimode Bose-Einstein condensates, subject 
to the action of a common resonant field, is analogous to a system of finite-level atoms in 
a common resonant electromagnetic field. A multitrap system can be formed, e.g., as an 
optical lattice with deep potential wells, incorporating many atoms around each lattice site. 
In the multitrap multimode condensate a high level of entanglement can be achieved. By 
varying the amplitude and frequency of the pumping resonant field, different regimes of 
evolutional entanglement can be realized. Moreover, by switching on and off the pumping 
field in various ways, it is feasible to create entanglement pulses of arbitrary length and 
composing arbitrary sequences of punctuated entanglement generation. Such a high level of 
admissible manipulation with and regulating of entanglement can, probably, be useful for 
information processing and quantum computing. 
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Appendix. Representation of coherent states in the 
Fock space. 



In this Appendix, we give, according to the referee's advice, an explanation of the main 
notations and definitions, used in the paper, which are related to the representation of 
coherent states in the Fock space. We mention here only the basic facts. Substantially more 
details, with an accurate mathematical foundation, cane be found in the very good book by 
Bcrczin [51]. Of course, there exist many other pubhcations expoimding this material, for 
instance, the classical textbook by Schweber [60]. Our notations are close to those employed 
in these books [51,60]. 

Recall, first, that the Fock space is a direct sum 



of n-particle Hilbert spaces 7i„ composed of symmetric or antisymmetric functions, according 
to whether the Bose-Einstein or Fermi-Dirac statistics are considered. For spinless bosons, 
treated in the paper, the wave function /„(ri, . . . , r„) is symmetric with respect to the 
permutations of the spatial variables rj, for all i = 1,2, ... ,n. The n-particle function 
/„(ri, . . . , Tn) pertains to the Hilbert space Tin- A state ip of the Fock space T is represented 
by the column 

/o 

/i(ri) 
/2(ri,r2) 
/3(ri,r2,r3) 



/n(ri,r2. 



It is clear that to write down each time such a column would be a too cumbersome way. 
Therefore one commonly shortens the notation, representing the state </? e ^ as 

(P^ [/„(ri,r2,...,r„)] , 

where only the generic element of the above column is written, with the index n = 0, 1, 2, . . . 
enumerating the rows of the column. Thus, the index n at the n-th row of the state </? tells 
that the n-particle wave function /„(ri, r2, . . . , r„) is implied. The scalar product of two 
states (fi, (f' of the Fock space T is defined as 



n=0 



where 



(/n, fn) = j fni^i, ^2, ■ ■ ■ , rn)/n(ri, r2, . . . , r„) dvi . . . dv^ 



is a scalar product in the n-particle Hilbert space 7i„. When a particular kind of wave 
functions /n(ri, r2, . . . , r„) is kept in mind, this is often emphasized by denoting the corre- 
sponding state in the Fock space as 



V=\f >^ [/n(ri,r2,...,r„)] . 
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Coherent states \rj > in the Fock space are defined as the eigenvectors of the annihilation 
field operator ip{r), so that 

ip{r)\ri > = ri{r)\ri > , 

where the coherent field ?7(r) plays the role of the eigenvalue. Being a member of the Fock 
space, a coherent state is also represented by a column 

Vo 

m(ri) 
?72(ri,r2) 



|77> = 



^n(ri,r2,...,r„) 



which is equivalently denoted by the accepted short-hand notation as 

k > = [?7n(ri,r2,...,r„)] . 

From the definition of the coherent states as the eigenvectors of the field operator and 
using the action of the latter yielding 

V'(r)|77 > = [VriTT77n(ri,r2, ...,r„,r)] , 

one gets the recursion relation 

Vn+1 77n+i(ri, r2, . . . , r„, r) = 77(r) ?7„(ri, r2, . . . , r„) . 

Solving this equation by repeated iterations, one finds 



77„(ri,r2, ...,r„J 



Vo 



The normalization condition < rjlrj >— 1 gives 

|77o| = expj- ^{r],r])^ . 
Assuming that the coherent field 77(r) is normalized to the number of particles N, one has 

(77, v) = J V*{''^)v{''^) dr — N . 
In this way, for a coherent state |?7 >, up to a global phase factor, one obtains 

1 



IT] > — e 



-N/2 



:^ V{ri)v{r2) 

:^ 77(ri)77(r2)77(r3) 

^ 77(ri)77(r2)...77(r„) 



This is exactly the coherent state (1) written down in the accepted short-hand notation. 

As is seen, to define the whole coherent state \r] >, one needs to find the coherent field 
r/(r). The latter, for the considered case of a dilute Bose gas with contact interactions, is a 
solution of the Gross-Pitaevskii equation. 
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Figure Captions 



Fig. 1. Regulated equidistant pulses of e2{t), formed by switching on and off the 
resonant field, with b = 0.7, so that B2{t) equals one during the time intervals At = 7.35 (in 
units of a~^), and it equals zero during the same intervals At — 7.35. Time is measured in 
units of a~^. 

Fig. 2. Nonequidistant pulses of e2{t), created by switching on and off the pumping 
field, with 6 = 0.7, at nonequal time intervals. Time is measured in units of a"^. 

Fig. 3. Regulated pulses of S2{t), for the same 6 = 0.7, as in Fig. 2, but for essentially 
different moments of switching on and off the pumping field. Time is measured in units of 
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FIG. 1. Regulated equidistant pulses of e2{t), formed by switching on and off the resonant field, 
with b = 0.7, so that £2it) equals one during the time intervals At = 7.35 (in units of a~^), and it 
equals zero during the same intervals At = 7.35. Time is measured in units of a~^. 
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FIG. 2. Nonequidistant pulses of £2{'t), created by switching on and off the pumping field, with 
b = 0.7, at nonequal time intervals. Time is measured in units of a~^. 
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FIG. 3. Regulated pulses of £2(^)1 for the same b = 0.7, as in Fig. 2, but for essentially different 
moments of switching on and off the pumping field. Time is measured in units of a~^. 
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